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Abstract. Error analysis of archived total 0 3 from total ozone mapping spectrometer 
(TOMS) (version 6) presented in earlier studies [Hudson and Kim, 1994; Hudson et al., 

1995] is extended to include scan angle effects. Daily total 0 3 maps for the tropics, from 
the period October 6-21, 1992, are derived from TOMS radiances following correction for 
these errors. These daily maps, averaged together, show a wavelike feature, which is 
observed in all latitude bands (10°N to 14°S), underlying sharp peaks which occur at 
different longitudes depending on the latitude. The wave pattern is used to derive both 
time-averaged stratospheric and tropospheric 0 3 fields. The nature of the wave pattern 
(stratospheric or tropospheric) cannot be determined with certainty due to missing data 
(no Pacific sondes, no lower stratospheric Stratospheric Aerosol and Gas Experiment 
(SAGE) ozone for 18 months after the Mt. Pinatubo eruption) and significant 
uncertainties in the corroborative satellite record in the lower stratosphere (solar 
backscattered ultraviolet (SBUV), microwave limb sounder (MLS)). However, the time- 
averaged tropospheric ozone field, based on the assumption that the wave feature is 
stratospheric, agrees within 10% with ultraviolet differential absorption laser Transport 
and Atmospheric Chemistry near the Equator-Atlantic) (TRACE A) 0 3 measurements 
from the DC-8 [Browell et al, this issue] and with ozonesonde measurements over 
Brazzaville, Congo (4°S, 15°E), Ascension Island (8°S, 15°W), and Natal, Brazil (5.5°S, 

35° W), for the period October 6-21, 1992. The derived background (nonpolluted) Indian 
Ocean tropospheric ozone amount, 26 Dobson units (DU), agrees with the cleanest 
African ozonesonde profiles for September-October 1992. The assumption of a totally 
tropospheric wave (flat stratosphere) gives 38 DU above the western Indian Ocean and 
15-40% disagreements with the sondes. Tropospheric column 0 3 is high from South 
America to Africa, owing to interaction of dynamics with biomass burning emissions 
[Thompson et al, this issue (a, b)]. Comparison with fire distributions from advanced very 
high resolution radiometer (AVHHR) during October 1992 suggests that tropospheric 0 3 
produced from biomass burning in South America and Africa dominates the 0 3 budget in 


spheric column 0 3 by subtracting stratospheric column 0 3 , 
obtained from the Stratospheric Aerosol and Gas Experiment 
(SAGE), from total ozone mapping spectrometer (TOMS) 
total 0 3 . Archived, gridded TOMS total 0 3 data with 1° X 
1.25° resolution (version 5, later version 6) was used for these 
studies. The tropospheric 0 3 residual maps showed enhanced 
0. 3 over the Atlantic off the coast of west Africa during the 
biomass burning season. This result drew great attention to 
tropical tropospheric chemistry. 

However, Thompson et al [1993], Hudson and Kim [1994], 
and Hudson et al [1995] have shown that the use of the ar- 
chived TOMS total 0 3 data from the current TOMS algorithm 
introduces an overestimation of total ozone over stratocumu- 
lus clouds by as much as 20 Dobson units (DU), and an un- 
derestimation of total ozone up to 15 DU when extremely high 
amounts of low-lying tropospheric O s are present. Since these 
conditions prevail over southern Africa and the adjacent At- 


the tropical southern hemisphere during the study period. 


1. Introduction 

Tropospheric ozone plays an important role in both tropo- 
spheric oxidation and global climate change and as a pollutant. 
Because ground-based measurements of tropospheric ozone 
give limited data in time and space, space-based measurements 
are in demand. Fishman and coworkers [Watson et al, 1991; 
Fishman, 1991; Fishman et al, 1991] developed the concept of 
a tropospheric “residual” to determine time-averaged tropo- 
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lantic during the biomass burning season, Hudson et al. [1995] 
suggested that the contrast in residual-derived tropospheric 
ozone between the Atlantic and the African continent is exag- 
gerated. In addition, unlike TOMS, which provides a daily 
global total ozone map with 50 X 50 km resolution, the spatial 
coverage of SAGE is limited to two latitude bands a day with 
15 measurements each. These are 25° longitude apart; latitude 
coverage is from 70° N to 70° S over a period of approximately 
1 month [McCormick et al., 1989]. Therefore significant differ- 
ences in spatial and temporal resolution between the two data 
sets can lead to an error in determining time-averaged tropo- 
spheric ozone for a short timescale. 

This paper presents a new method to overcome these diffi- 
culties in determining time-averaged tropospheric column 0 3 
in the tropics using only TOMS radiances. Since the archived 
TOMS total 0 3 (version 6) has errors associated with clouds, 
low retrieval efficiency of tropospheric ozone, Pinatubo aero- 
sols [Bhartia et al., 1993], and the determination of the space- 
craft attitude [Aerospace Engineering and Research Associates 
( AERA ), 1990], we have used TOMS radiances to correct these 
errors (section 2). Section 3 shows a new method to derive 
time-averaged stratospheric and tropospheric 0 3 fields over 
the entire tropical band using only TOMS data and sondes. 
The error analysis of this method is discussed in section 4. The 
time-averaged tropospheric column 0 3 derived using the new 
method is then compared with ozonesonde and ultraviolet 
differential absorption laser (uv-DIAL) ozone measurements 
over the South Atlantic basin (section 5). Finally, features of 
the tropospheric 0 3 budget over the tropics during October 
1992 are discussed (section 6). 

The period of study is October 6 to 21, 1992, which is of 
interest because of South Atlantic field experiments: Southern 
Tropical Atlantic Regional Experiment (STARE) Southern 
African Fire Atmospheric Research Initiative (SAFARI) and 
Transport and Atmospheric Chemistry near the Equator- 
Atlantic (TRACE A), in September and October 1992 [An- 
dreae et al., 1994a; Fishman et al., this issue (a)]. 

2. Corrections to Archived Total Ozone 
From TOMS 

The TOMS instrument scans across a track perpendicular to 
the spacecraft velocity vector in 35 sequential scan positions 
[McPeters et al., 1993]. Scan positions 1, 18, and 35 correspond 
to the. highest scan angle to the right, the nadir angle, and the 
highest scan angle to the left, respectively. The coverage for 
one cross track corresponds to about 26° longitude at the 
equator, and thus it takes about 14 orbits per day to cover the 
entire globe. The 14 orbits actually cover about 365°, so the 
nadir track shifts 5° westward each day. Figure la shows ar- 
chived total 0 3 from the high-density tapes over the tropics for 
October 1 and 2, 1992. The total 0 3 distribution shows three 
maxima at 178°E, 156°W, and 132°W, about 26° apart, which 
shift westward about 5° on the following day. Because the 
number of maxima, the interval between maxima, and the 
motion of maxima exactly match, respectively, the number of 
the satellite orbits, TOMS coverage for each orbit, and orbit 
shift, these maxima most likely represent an apparent TOMS 
scan bias. 

Two factors contribute to this bias: (1) an error due to 
scattering by stratospheric aerosols in the stratosphere arising 
from the Mt. Pinatubo eruption [Bhartia et al., 1993] and (2) an 
error in the Nimbus 7 satellite roll angle [AERA, 1990]. As the 


satellite roll angle changes, so does the effective scan angle at 
the Earth. The error in the derived total 0 3 is more at the 
extreme off-nadir positions (1 and 35) than at the nadir posi- 
tion (18) because path length changes as the secant of the 
angle [McPeters et al., 1993]. In order to determine the error in 
the derived total 0 3 due to the attitude error, we averaged the 
measured total 0 3 under clear sky conditions (reflectivity less 
than 0.2) as a function of scan position over the tropics in 
October 1989 prior to the Mt. Pinatubo eruption. Figure 2a 
shows the averaged total 0 3 in each scan position relative to 
that in scan position 18, where the 0 3 error is negligible [Mc- 
Peters et al., 1993]. Total 0 3 is underestimated at the right scan 
position (1-17) and overestimated at the left scan position 
(19-35). 

The error due to stratospheric aerosols arises because the 
TOMS algorithm assumes that only Rayleigh scattering occurs. 
In order to account for this error, we averaged the measured 
total 0 3 for clear sky conditions as a function of scan position 
from October 6 to 21, 1992. Because the averaged total 0 3 in 
1992 contains both the attitude error and the aerosol error, the 
attitude error derived above is first subtracted from the aver- 
aged total 0 3 . The error due to aerosols in the stratosphere has 
been shown to be negligible at scan position 1 [Bhartia et al., 
1993]. Hence we can derive the aerosol error by subtracting the 
corrected total 0 3 at scan position 1 from the total 0 3 at the 
other scan positions (Figure 2b). Figure 2c is the total error as 
a function of scan angle. The total error is subtracted from the 
derived total 0 3 to obtain corrected total 0 3 (see Figure lb). 
Note that the scan bias has been eliminated. 

3. A New Algorithm to Determine Time- 
Averaged Stratospheric and Tropospheric 
Column Ozone 

Ideally, in order to retrieve the tropospheric column 0 3 
using only TOMS radiances, we require a knowledge of the 
stratospheric column 0 3 , the tropospheric ozone profile, and 
stratospheric ozone profile. However, Hudson et al. [1995] and 
Kim [1995] have shown that in the tropics the algorithm does 
not depend on the stratospheric ozone profile. The following 
procedures have been used to derive corrected total 0 3 , time- 
averaged stratospheric column 0 3 , and tropospheric column 
0 3 using measured tropospheric 0 3 profiles. 

3.1. Derivation of Total 0 3 , Xl tot From TOMS Radiances 

The current TOMS algorithm in the tropics uses a series of 
lookup tables based on standard profiles of 225, 275, and 325 
DU. The amount assumed to be in the tropospheric is approx- 
imately 10% of the total amount, thus implying a coupling of 
stratospheric and tropospheric ozone. Table 1 shows climato- 
logical tropospheric column 0 3 for tropical and subtropical 
sites, including seven that were operational during SAFARI 
and TRACE A [Baldy et al., this issue; Kirchhoff et al., this 
issue; Thompson et al., this issue (a, b)]. For non-SAFARI and 
TRACE A sites, September to November means are 26 DU at 
Hilo (19°N, 155°W) and Samoa (14°S, 170°W) and 24 DU at 
Manaus (3°S, 60°W) [Komhyr et al., 1989; Watson et al., 1991]. 
Note that these values are similar to the minima observed 
during SAFARI/TRACE A at sites that recorded clean back- 
ground soundings: Brazzaville, Cuiaba, and Okaukuejo. Val- 
ues of —20-30 DU are also typical for the tropical South 
Atlantic basin sites during the non-biomass burning season 
(Table 1). Weller et al. [1996] show mid-Atlantic ozone profiles 
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(a) Archived TOMS Total Ozone (DU) on Oct 1 , 1 992 




Derived TOMS Total Ozone (DU) on Oct 2, 1992 



Figure 1 . (a) Archived total 0 3 for the tropics from high-density TOMS tape, October 1 and 2, 1992. (b) 

Total 0 3 after correction for error due to a scan bias. Five DU intervals shown. This Pacific region was 
selected because tropospheric ozone is a minimum and it introduces no scan angle artifact into total ozone. 
Thus features due to attitude and Mt. Pinatubo aerosol effects are most pronounced. 


from shipboard launches during October and November 1993. 
Background tropospheric ozone column (0°-10°N, tropopause 
at 16-17 km) is 27 DU. In general, the SAFARI/TRACE A 
sondes show no evidence that tropospheric and stratospheric 
total ozone amounts are consistently coupled, in agreement 
with the shipboard record [Weller et aL, 1996]. 


In order to obtain total ozone amounts from the measured 
radiances we have followed the same lookup table procedure 
as used in the TOMS algorithm, but new tables were generated 
from a set of profiles with varying stratospheric total ozone 
amounts and a fixed tropospheric column 0 3 of 26 DU. A 
description of the algorithm is given by Hudson et al. [1995]. 
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Ozone Error Due to Aerosol, Satellite Attitude, and Tropospheric Ozone Profile 



10 20 30 


Satellite Attitude, and 
Tropospheric Ozone Profile 


Pinatubo Aerosol 


Total Error 


Scan Position (1 - 35) 

Figure 2. Ozone error due to (a) the error in scan angle and (b) Mt. Pinatubo aerosols, (c) Total error due 
to the error in scan angle and Mt. Pinatubo aerosols. 


Since the retrieval efficiency for tropospheric 0 3 is less than 
that for stratospheric 0 3 [Hudson et ah, 1995], the algorithm 
underestimates the total 0 3 when the actual tropospheric col- 
umn 0 3 is greater than that assumed in the tables. The cor- 
rection for this underestimation is discussed in section 3.2. The 
derived total 0 3 over cloud-free regions based on the 26-DU 
tropospheric 0 3 assumption, and averaged over 2° longitude 
and 1° latitude bands for 16 days in October 1992, is shown in 
Plate 1. 

3.2. Derivation of Time-Averaged Stratospheric and 
Tropospheric O a 

Figure 3a shows the longitudinal variations of total 0 3 for 
the latitude band from 2°S to 2°N. Figure 3b shows a plot of the 
longitudinal variation for latitude bands at 14°S, 8°S, and 10°N, 
normalized to that shown in Figure 3a. The longitudinal vari- 
ation of total 0 3 at the equator is selected because the varia- 
tion of stratospheric ozone with longitude has been measured 


to be smooth at the equator [Fishman et al., 1990] and is 
representative of the variation in other latitude bands. The 
areas from 180°W to 90°W (over the Pacific) and from 50°E to 
100°E (over the Indian) are selected for this normalization 
procedure because they are least affected by anthropogenic 
activity [Fishman and Larsen, 1987; Fishman et al., 1990; Kom- 
hyr et al., 1989]. Note that the shapes of the total 0 3 fields track 
each other in the latitude bands from 14°S to 1Q°N, except from 
Africa to South America. The composite plot shown in Figure 
3b shows a wavelike feature, which is observed in all the lati- 
tude bands, underlying sharp peaks which occur at different 
longitudes depending on the latitude. The wave feature is quite 
distinct from 10°N to 14°S. Beyond these latitudes the wave 
shows a more complicated structure. 

We have performed a similar analysis for the first half of 
each month during 1992 [Hudson et al., 1996] and for early 
October 1989. In the first case, all months in 1992 show a wave 
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Table i. Climatological Integrated Tropospheric Ozone for Non-Burning Season 
(March- April-May) and Minimum Value During SAFARI/TRACE A for South Atlantic, 


Brazil, African, and Indian Ocean Stations 


Station 

Location 

March-April-May 
Mean, DU 

Sept.-Oct. 1992 
Minimum, DU 

Ascension Island 

(8°S, 15°W) 

31 

37 

Brazzaville 

(4°S, 15°E) 

33 

24 

Goiania, Brazil 

(16°S, 49°W) 

22 

n/a 

Natal, Brazil 

(5°S, 35°W) 

27 

40 

Cuiaba, Brazil 

(15.5°S, 56°W) 

n/a 

28 

Okaukuejo, Namibia 

(19°S, 16°W) 

n/a 

23 

La Reunion 

(22°S, 56°E) 

28 

30 

Irene, South Africa 

(26°S, 28°E) 

29 

29 

Tutuila, American 

(14°S, 170°W) 

21 

n/a 

Samoa 
Hilo, Hawaii 

(20°S, 155°W) 

41 

n/a 


Ozone is in Dobson units to 16 km. 


structure similar to Figures 3a and 3b, from 4°N to 5°S; beyond 
this band the definition of wavelike structures depends on the 
extent of the intrusion of midlatitude waves into the tropics. 
For October 1989 the wave structure is also strongly evident 
(Figure 3c), despite the fact that total ozone in the 10°N-to- 
14°S band is 4 DU less than its October 1992 counterpart. 

For nearly all latitude bands shown in Figures 3b and 3c 
elevated total ozone occurs between 60°W and 40°E. North of 
the equator these correspond to regions of biogenic and in- 
dustrial emissions conducive to ozone formation. The ozone 
precursor source strengths are not known, but if the enhance- 
ment above the normalized background is entirely due to tro- 
pospheric ozone, the amounts seem high. (We have not been 


able to obtain in situ ozone data north of the equator for any 
of the periods studied, but we derive tropospheric ozone col- 
umns in excellent agreement with sondes south of the equa- 
tor). Upper tropospheric ozone formation, perhaps lightning- 
NO-induced, is an alternative explanation. For October 1992, 
forward trajectories corresponding to altitudes of convective 
tower outflow (9-12 km), could carry air masses elevated in 
ozone or ozone precursors from Africa toward the Gulf of 
Guinea [Thompson et al., this issue (a)]. South of the equator, 
aircraft sampling on TRACE A showed that elevated ozone 
and ozone precursors were abundant over Africa and the east- 
ern and central Atlantic [Gregory et al., this issue; Browell et al., 
this issue]. 



Plate 1. Derived time-averaged total column 0 3 over the tropics. Color bar represents column 0 3 in 5-DU 
intervals. 
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Time-averaged Total 0 3 Averaged Between 2°S and 2°N 
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Normalized Time-averaged Total 0 3 (6-21 Oct 1992) 
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Figure 3. (a) Time-averaged total 0 3 averaged over latitude bands from 2°S to 2°N for October 6-21, 1992. 

(b) Normalized time-averaged total 0 3 and background total 0 3 (solid line) over the tropics. Each grid point 
is averaged over 2° longitude and 1° latitude, (c) Normalized time-averaged ozone and background total O s 
for October 6-21, 1989. 
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Figure 3. (continued) 
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The wave structure in Figure 3 suggests an approach to 
deriving tropospheric ozone. Consider that total ozone, fl tot , 
consists of a highly variable amount above the normalization 
curve (Figure 3), plus the amount indicated by the latter, the 
background amount of ozone, fl bck . Then the difference be- 
tween those two quantities, corrected with the retrieval effi- 
ciency for tropospheric ozone (e, less than unity for high con- 
centrations of ozone below 500 mbar), is 

a tr = [tttot T" O tr> bck (1) 

where fl tr bck is a background amount of tropospheric ozone. 
All quantities are in Dobson units. Now what is H tr bck ? The 
minimum amount, as indicated above from the soundings, is 26 
DU independent of longitude. In that case the wave structure 
is assumed to be in the stratosphere and fl bck = D strat + 26 
DU; fi bck at a given latitude is actually a modified parameter, 
derived by renormalizing the originally derived fl bck to the 
wave 1 curve in Figure 3b. This assumption leads to Plate 2 for 
the time-averaged tropospheric ozone field for October 6-21, 
1992. Since the retrieval efficiency depends on the tropospheric 
0 3 profile, we use averaged profiles derived from Ascension 
Island, Brazzaville, and Okaukuejo ozone soundings during 
September and October 1992 (Table 1). Calculated time- 
averaged retrieval efficiency (e) is about 0.75 under clear sky 
conditions and about 1.0 over marine stratocumulus clouds in 
the tropics. 

An alternative is to assume a longitudinally invariant strato- 
spheric ozone amount, which is not very different from what 
SAGE usually shows, though SAGE was not available for late 
1992, owing to the Mt. Pinatubo eruption. This was the ap- 
proach used by Hudson et al. [1995], and it effectively ascribes 


the wave structure in Figure 3 to tropospheric ozone. Because 
SAGE data were not available, we have assumed that the 
stratospheric ozone field can be obtained by subtracting 26 DU 
from the total ozone near 180° longitude. Thus (1) is still used 
to derive H tr , but O bck = D tr bck plus a fixed amount assumed 
for fl strat . Plate 3 shows the averaged tropospheric ozone field 
for October 6-21, 1992, with this assumption. 

3.3. Equatorial Wave Structure 

A number of studies have considered wave structures in 
equatorial 0 3 that are apparent when total ozone and the 
deviation from zonal mean total ozone are viewed as a function 
of longitude. A Kelvin wave [Andrews et al., 1987; Tsuda et al., 
1994], with wave number 1 or 2 and with both 7-day and 15- to 
20-day periods, has a 2- to 4-DU magnitude [Ziemke and 
Stanford, 1994]. Another contribution to total ozone is a mixed 
Rossby-gravity wave with wave number 4, and a 4- to 5-day 
period, also of magnitude 2-4 DU [Stanford and Ziemke, 
1993]. Shiotani and Hasebe [1994] found no evidence of a wave 
in SAGE ozone corresponding to that seen in TOMS and 
concluded that the wave must be primarily tropospheric. How- 
ever, SAGE ozone profiles are >30% uncertain below 100 
mbar, and a wave pattern would be difficult to isolate. Table 2 
summarizes the SAGE-derived ozone thickness in 2-km layers 
in the upper troposphere-lower stratosphere, for all observa- 
tions between 10°N and 10°S in 1989-1990. The second col- 
umn of Table 2 is ozone integrated from top of atmosphere 
down to the specified altitude. 

Comparison of TOMS total ozone and stratospheric ozone 
from the microwave limb sounder (MLS) on UARS) also sug- 
gests a primarily tropospheric wave [ Ziemke et al., 1996], al- 
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time-averaged tropospheric column 03 (wave-one) 



22 25 27 50 55 56 59 4-2 45 4^ 52^ ^ 55 58 61 64 67 70 

Column 03 in DU 



Plate 2. Derived time-averaged tropospheric column 0 3 over the tropics, assuming that the wavelike 
feature in Plate 1 is in the stratosphere. Color bar represents column 0 3 in 3-DU intervals. 


though ozone is > 40 % uncertain below 68 mbar, (E. Fishbein, sondes (September-October for 1986-1989) implies that the 
Jet Propulsion Laboratory, UARS Science Team documenta- wave is largely tropospheric, but no Pacific sondes were avail- 
tion, 1994). A comparison of TRACE A Ascension sondes able for an actual differencing of Atlantic-Pacific ozone pro- 
(September-October 1992) with a climatology of Samoa files in 1992. Preliminary analysis with version 7 TOMS ozone 


time-averaged tropospheric column 03 (flat stratospheric 03) 



22 25 27 30 33 36 39 42 45 48 52 55 58 61 64 67 70 

Column 03 in DU 



Plate 3. Same as Plate 2, except assuming a flat stratosphere, that is, a tropospheric wave. The latter is 
similar to a “residual method.” Color bar represents column 0 3 in 3-DU intervals. 
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including all months of 1992 shows that the wave pattern is 
present year-round (discussed by Ziemke et al. [1996] and C. N. 
Hsu et al., The effect of an improved cloud climatology on the 
TOMS total ozone retrieval, submitted to Journal of Geophysi- 
cal Research, 1996] and that contributions of stratospheric and 
tropospheric ozone to the wave pattern vary from month to 
month [Hudson et al., 1996]. In section 5 it is shown that 
assuming a stratospheric wave gives better agreement with 
TRACE A sondes over Brazil, Ascension, and Africa in Oc- 
tober 1992. 

4. Error Analysis of the New Method 

By far the largest source of error is in the choice of the 
stratospheric ozone field. If the wave 1 pattern is in the strato- 
sphere, then an error arises from uncertainties in drawing the 
best-fit underlying pattern in Figure 3b. The shape of this curve 
is based largely on the longitudinal variation of total 0 3 at the 
equator. The magnitude of this error can be estimated from 
the standard deviation of the total 0 3 at a given latitude from 
the best-fit curve. The standard deviations, calculated from 
180°E to 80°W and from 60°E to 100°E, which are the regions 
with less anthropogenic activity, are less than 2 DU. This leads 
to an error in tropospheric 0 3 of 3 DU. 

The second error arises from uncertainty in the retrieval 
efficiency. Hudson et al. [1995] suggested that this error is less 
than 2 DU in derived tropospheric column 0 3 on a daily basis. 
This error occurs due to the difference between actual and 
assumed tropospheric 0 3 profile shape on a daily basis. Pre- 
sumably, this error is less over the period of 15 days because we 
used averaged 0 3 profiles measured from ozone soundings 
during the same period as our case study. 

The third error will occur where background tropospheric 
0 3 differs from the assumed background 0 3 of 26 DU. The 
difference would propagate directly to an error in derived 
tropospheric 0 3 . For example, if the actual background tropo- 
spheric 0 3 is 31 DU instead of 26 DU, the derived tropo- 
spheric 0 3 will be underestimated by 5/e, which are about 7 
DU under clear sky and 5 DU under cloudy sky conditions, 
respectively. Looking at the minimum 0 3 values during SA- 
FARI/TRACE A for Brazzaville, Cuiaba, and Okaukuejo (Ta- 
ble 1), it seems unlikely that this error could be larger than 5 
DU. Indeed, we have tried to use a moderate value of “back- 
ground” ozone to avoid exaggerating regional variations. 

5. Intercomparison of TOMS-Derived 
Tropospheric 0 3 With SAFARI/TRACE A Data 

Table 3 shows a comparison between total tropospheric 
ozone derived from the sonde measurements and that derived 


Table 2. SAGE-Derived Ozone Thicknesses in 2-km 
Layers in the Upper Troposphere-Lower Stratosphere 


Altitude, 

km 

Ozone Column 
to Altitude, DU 

Incremental 
Ozone, DU 

lo- Uncertainty in 
Incremental 
Ozone, DU 

19 

230.35 (619) 



17 

236.73 (617) 

15.4 

2.63 

15 

241.94 (581) 

6.39 

2.43 

13 

244.90 (411) 

2.96 

7.20 


Thicknesses are for all observations between 1°N and 1°S in 1989- 
1990. The second column is ozone integrated from top of atmosphere 
down to the specified altitude. Numbers in parentheses are numbers of 
observations. DU, Dobson unit. 


Table 3. Sonde and TOMS-Derived Tropospheric Ozone 
From October 6 to 21, 1992 

Ozone, DU 


Location (latitude, longitude) 

Sonde- 

Derived 

TOMS- 
Derived 
(Plate 2) 

TOMS- 
Derived 
(Plate 3) 

Brazzaville (40°S, 150°E) 

40(7) 

38 

57 

Ascension Island (8°S, 15°W) 

52(6) 

48 

58 

Natal (5°S, 35°W) 

48 (4) 

46 

62 

Cuiaba (15.5°S, 56°W) 

41 (3) 

40* 

50* 

Samoa (14°S, 170°W) 

26 (climatology) 

25 

28 


“Plate 2” refers to identifying the wave as stratospheric; “plate 3” 
refers to identifying the wave as tropospheric. Numbers in parentheses 
indicate the number of ozonesonde measurements. 

* Comparisons from southern limit of Plates 2 and 3: 14°S, 56°W. 


from the new method for the two different assumptions about 
the nature of the wave 1 pattern. The number in the paren- 
thesis in the second column is the number of ozonesonde 
measurements for the period October 6-21, 1992. Since deri- 
vation of ozone from this method is limited to 14°S-10°N, the 
ozone value at 14°S, 56°W is given for comparison with Cuiaba, 
Brazil. The value at Samoa is the climatological tropospheric 
ozone from September to October [Komhyr et al., 1989]. Best 
agreement between the tropospheric ozone derived from the 
new method and the sondes is for the assumption that the wave 
one pattern is in the stratosphere. 

Figure 4 shows a composite of integrated tropospheric col- 
umn ozone taken by the Langley uv-DIAL instrument on 
board the NASA DC-8 during the TRACE A aircraft mission 
in October 1992 [Browell et al., this issue]. The uv-DIAL mea- 
sured ozone profiles above and below the DC-8 simulta- 
neously, and the integrated column amount was obtained by 
combining above- and below-aircraft ozone amounts. Flights 
13-18 were west of Africa over the Atlantic, with flights 13-15 
originating and returning to Windhoek, Namibia, between Oc- 
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Figure 4. Uv-DIAL integrated tropospheric ozone from 
flights over southern Africa and the tropical Atlantic; flights 
occurred from October 6 to 22, 1992. Flight tracks are given by 
Fishman et al. [this issue (a)]. Longitudinal range of data dis- 
played is 20°W to 30°E. Sun angle prevented collection of data 
between 4°S and 10°N [from Browell et al., this issue]. 
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Table 4. Area- Averaged Tropospheric 0 3 and Budget 
Over the Tropics 


Tropospheric 0 3 
Column 


Location 

Longitude, Latitude 

2, DU 

Budget, Tg : 

South America 

75°-40°W, 0°-14°S 

39 

1.6 

Atlantic 

40°W-12 q E, 0°-14°S 

45 

3.5 

Southern Africa 

12°-50°E, 0°-14°S 

37 

1.5 

Area 1 

12 o -50°E, 0°-5°S 

31 

0.3 

Area 2 

12°-50°E, 5°-10°S 

38 

0.6 

Area 3 

12°-50°E, 10°-14°S 

44 

0.6 


Areas within indicated latitude-longitude ranges: South America, 
6.1 x 10 6 km 2 ; Atlantic, 9.0 X 10 6 km 2 ; Africa, 6.6 X 10 6 km 2 . 

*1 tg is equal to 1 X 10 12 g. 


tober 14 and 18, 1992. Flight tracks are given by Fishman et al. 
[this issue (a)]; the longitudinal range of the data displayed in 
Figure 4 is 20°W to 30°E. The mean between 4°S and 15°S is 50 
DU, in good agreement with the 45-DU amount from the 
TOMS-derived tropospheric ozone method (Table 4). On the 
one continental flight (Flight 10, October 6, 1992) that reached 
to biomass burning regions of Africa (northern Zambia, 10°S, 
30°E), column integrated ozone was 51-53 DU. This is just 
within the range of our map in Plate 2, where the ozone 
column is 51-54 DU. The map shown in Plate 3 does not agree 
with the uv-DIAL results. Values of 50-55 DU over the east- 
ern Pacific and about 70 DU over south-central Africa in Plate 
3 are too high. The TOMS/soIar backscattered ultraviolet 
(SBUV) method [Fishman et al., this issue (b)] of deriving 
tropospheric ozone also does not reproduce the northern 
Zambian 0 3 maximum seen in the uv-DIAL. 

In summary, both the sondes and the uv-DIAL measure- 
ments support the tropospheric ozone map given in Plate 2. 
The remainder of our analysis is based on these data. 

6. Analysis of Derived Tropospheric Column 
Ozone 

The overall features from Plate 2 are (1) very high ozone 
from South America to Africa, (2) relatively low ozone off the 
western coast of South America and from the Indian to the 
west Pacific Ocean, (3) relatively enhanced ozone over north- 
ern Australia and New Guinea, and (4) a plumelike structure 
extending from the eastern Pacific to central Pacific. We dis- 
cuss the Atlantic and Pacific panels of Plate 2 separately. 

6.1. Tropospheric 0 3 From Eastern South America 
Through Africa 

Tropospheric 0 3 of about 40 DU is observed over northwest 
central Africa near the Intertropical Convergence Zone 
(ITCZ), (equator to 10°N). The source of elevated tropo- 
spheric 0 3 over this region is uncertain. Because the biomass 
burning season in northern Africa is from December to Feb- 
ruary [Marenco et al., 1990; Cros et al., 1991], local burning over 
the savanna (8°N to 15°N) is probably not the cause of en- 
hanced tropospheric 0 3 . Biogenic sources of ozone precursors 
(e.g., hydrocarbons, CO, and NO) [Zepp et al., this issue; Harris 
et al., this issue] from the rain forest (8°N to 3°S and 15°W to 
30°E) could contribute 2- to 3-DU boundary layer ozone. 
Ozone from lightning NO near the ITCZ and/or ozone formed 
following convective outflow was observed during TRACE A 


Brazilian sampling [Pickering et al., this issue]. In general, Ol- 
son et al. [this issue] and Thompson et al. [this issue (a), 1996] 
note that in September and October 1992, the highest tropo- 
spheric ozone recorded at Natal and Ascension occurred when 
upper level westerlies transport air masses from convectively 
active regions during the latter part of the South American 
biomass burning season. However, back trajectories run with 
the Goddard isentropic trajectory model from —300 and 700 
mbar show that the 0 3 feature of 50 DU near 3°N, 8°W orig- 
inates from northwest central Africa. There is no trace gas data 
available for this region, but industrial pollution could be 
present [Sieber et al., 1995; Thompson et al., 1996, this issue 
(a)]. Northeast central Africa, in contrast, is under the influ- 
ence of clean air transported from the Indian Ocean (see mean 
winds in Figure 5) and does not reflect significant burning 
input because most African savanna burning during October 
1992 occurs between 10° and 20°S [Justice et al., this issue]. 

Figure 6 shows southern African fire counts for October 
1-15, 1992, in southern Africa. Southeast central Africa is close 
to the intensive northern Zambian burning region, with a tro- 
pospheric 0 3 column up to 60 DU. This suggests that 0 3 is 
rapidly formed from fresh precursors in the burning regions. 
Model calculations using SAFARI and TRACE A aircraft 
measurements of 0 3 and O s precursors (CO, NO, nonmethane 
hydrocarbons) in the lower troposphere near active burning 
show that a typical cross section of air influenced by biomass 
burning (from boundary layer to tropopause) produces 10-15 
ppbv 0 3 /day (—5 DU column integrated). Furthermore, except 
in the boundary layer, the 0 3 photochemical lifetime is several 
weeks. This lifetime, coupled with lifetimes of stable, recircu- 
lating layers of or greater than 1 week, over southern Africa 
[Garstang et al., this issue; Tyson et al., 1996] pauses ozone to 
accumulate in this region. The persistent anticyclone that 
causes this recirculation means that south of 18°S, only 4% of 
air masses exits Africa toward the Atlantic Ocean [Garstang et 
al., this issue]. 

In the upper troposphere over Africa and the Atlantic, 
TRACE A measurements suggest that there is net 0 3 produc- 
tion rate at 1-5 ppbv 0 3 /day. This corresponds to tropospheric 
0 3 column production of 1-2 DU/day, a feature that shows up 
in time versus longitude (Hovmoller) plots of TOMS total 0 3 
[Thompson et al., this issue (b), 1996]. Thus a high tropospheric 
0 3 column over the eastern Atlantic reflects transport of low- 
level to midlevel ozone from Africa, supplemented by a steady 
rate of upper tropospheric ozone formation during transit 
from both continents. This is seen in the high south Atlantic 
ozone shown in Plate 2 and the wind fields shown in Figure 5. 

Biomass burning-linked 0 3 coincident with middle tropo- 
spheric flows from both South America and Africa [Andreae et 
al., 1994b] were also observed on western Atlantic flights dur- 
ing Chemical Instruments and Test Evaluation-3 (1989). 
Thompson et al. [this issue (a)] posed the following question: 
How much of the tropospheric ozone above background level 
(“excess ozone”) over the south Atlantic basin during SAFARI/ 
TRACE A was due to African origins and how much was due 
to South American origins? By running back trajectories from 
ozone-enriched layers observed on every Natal, Ascension, and 
Brazzaville sounding, it was determined that 20-25% of excess 
ozone at Natal and Ascension was from South America, with 
70-75% from Africa. Up to 10% of the excess had recirculated 
so long over the Atlantic that the ultimate continental origin 
was indeterminate. All the Brazzaville soundings corresponded 
to African origins. 
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Total Fire Counts, Oct 1-15, 1992 
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Figure 6. Satellite-derived fire counts for October 1-15, 1992 [Justice et al., this issue]. 


The fact that both South America and Africa apparently 
supply the Atlantic with O s and 0 3 precursors explains the 
high 0 3 column over the ocean. The ozone amounts over 
South America, Africa, and the Atlantic can be used to esti- 
mate a regional budget for 0 3 from biomass burning. Table 4 
gives the area for each region. We designate as 2 (additional 
averaged tropospheric ozone) the tropospheric 0 3 beyond back- 
ground level and convert to teragrams (10 12 grams) of ozone for 
each area (area in square centimeters) as follows: 

(2 ~ 26) X 2.69 X 10 16 X area X 48/6.023 X 10 23 

where 26 DU is background tropospheric column 0 3 (1 DU = 
2.69 X 10 16 (molecules/cm 2 )), 48 is the molecular weight of 0 3 
in grams per mole, and 6.023 X 10 23 is Avogadro’s number, in 
molecules per mole. This is the accumulated 0 3 budget above 
background which is a function of the rate of photochemical 
0 3 production from biomass burning and biogenic precursors 
as well as the lifetime of 0 3 . 

We calculate a total ozone burden for southern Africa for 
October 6-21, 1992. Since this study does not extend south of 
14°S, we must estimate tropospheric 0 3 south of 14°S to find 
the total budget produced by biomass burning. From remote 
sensing data over Africa [Justice et al., this issue] fire counts 
north of 14°S are twice what they are south of 14°S. It is 
reasonable to assume that the same ratio holds in South Amer- 
ica (C. Justice, personal communication, 1994). For a band 
from 0° to 30°S between Africa and South America, this results 
in a total tropospheric 0 3 burden for this region of 9.9 Tg 
(Table 4). If an average 0 3 formation rate of 2-3 DU/day 
[Thompson et al., this issue (b)] is assumed, the 0 3 formation 
rate from South America, across the Atlantic, to Africa is 1.4 
—2.1 Tg/day, with about a 6-day lifetime (9.9(Tg)/{1.4 — 2.1 


(Tg/day)}). Therefore the burden of tropospheric 0 3 produced 
from biomass burning in October is 51 Tg (9.9 (Tg) X 1 
(month)/6 (days)). 

In order to derive a total burden from southern tropical 
African and South American biomass burning, we have to 
calculate the burden in each month (June to September) rel- 
ative to that in October. For this calculation we use the records 
of integrated tropospheric 0 3 in Brazzaville, Ascension Island, 
and Natal [Cros et al., 1992; Fishman et al., 1992]. These show 
tropospheric 0 3 beyond background level to be approximately 
the same in August through October, with October about 1.5 
times greater in tropospheric 0 3 than in June and July. Thus 
we compute for the total burden of biomass-burning-derived 
tropospheric 0 3 (June to October) from South America, the 
Atlantic, and southern Africa: 

51X3 + 2/3 

X {9.9 (Tg) X (2 months)/(6 days)} ~ 220 Tg 0 3 

Because surface measurements of tropospheric 0 3 are lower 
and the period of biomass burning season in northern Africa 
(December to February) is shorter than those in southern 
Africa [Cros et al, 1991], the total burden of tropospheric 0 3 
in northern Africa is probably less than in southern Africa. 

Andreae [1991] suggests that the 0 3 burden from biomass 
burning over South America and Africa is 230 Tg 0 3 /yr with an 
error of factor of 2. Lelieveld et al. [1996], in a global (three- 
dimensional) model calculation of 0 3 sources and sinks, com- 
putes a net 0 3 loss in the tropics due to high UV insolation and 
humidity: via the H 2 0 + 0( 3 £>) reaction, 0 3 is destroyed. 
Biomass burning inputs do, however, result in formation of 144 
Tg 0 3 /yr. 
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6.2. Tropospheric 0 3 Over Western South America 
and the Pacific 

Plate 2 shows that averaged tropospheric column 0 3 (28 
DU) over the east Pacific Ocean is slightly higher than that (26 
DU) over the central Pacific Ocean, but more in situ observa- 
tions are required to say this with certainty. If there is elevated 
tropospheric 0 3 over the eastern Pacific due to transport from 
South America, it would not be as efficient as the transport 
from Africa to the Atlantic because of the high barrier of the 
Andes. These mountains also contribute to the very low 0 3 
observed over western South America. The Pacific Exploratory 
Mission-Tropics (PEM-Tropics) aircraft mission, scheduled 
for August-October 1996, will build a tropospheric ozone cli- 
matology for the tropical Pacific, with one or two ozone sound- 
ing sites being added to Samoa. 

A tropospheric 0 3 column as high as 40 DU is observed near 
140°E, where northern Australia and New Guinea are located. 
This region is affected by biomass burning from these nations 
and from Indonesia and Borneo [Hurst et al., 1994]. The in- 
tensities in Plate 2 suggest that this area contributes much less 
to the tropical 0 3 budget than Africa or South America. The 
lowest tropospheric column 0 3 is observed from the Indian to 
the far western Pacific Ocean. 

7. Conclusion 

This paper continues a series of studies aimed at deriving the 
tropospheric 0 3 column in the tropics from high-density TOMS 
radiances. Errors inherent in archived TOMS 0 3 data have been 
corrected: cloud effects (cf analysis of Hudson and Kim [1994] 
with archived TOMS of Thompson et al. [1993]), inefficiency of 
lower tropospheric 0 3 retrieval [Hudson et al., 1995], and sat- 
ellite attitude error and aerosols (this work). The present error 
analysis emphasizes the need to use the TOMS radiances in 
the high-density tapes to derive tropospheric column 0 3 . 

A new method has been developed to derive the time- 
averaged stratospheric ozone field over the tropics, a necessary 
precursor to determining the tropospheric column ozone. 
From examination of the total ozone field, an underlying wave 
1 pattern, which could be ascribed to a stratospheric or tropo- 
spheric wave, has been identified. For the period of October 
6-21, 1992 (during the TRACE A experiment), the tropo- 
spheric ozone map based on a stratospheric wave assumption 
gives excellent agreement (within 10%) with the DC-8 aircraft 
ozone data and with tropical ozone soundings at Brazzaville, 
Ascension Island, and Natal (Brazil). Further elucidation of 
the wave requires 0 3 profile measurements from surface to 
lower stratosphere. Our method is valid only in the tropics 
where stratospheric ozone is relatively invariant over several- 
week periods of time. The fact that the sondes are 2-4 DU 
greater than our map suggests that we have overestimated 
stratospheric ozone by this amount, but it is within the uncer- 
tainty of the method. Beyond the tropics, derivation of accu- 
rate tropospheric 0 3 from TOMS requires detailed knowledge 
of the stratospheric ozone profile. 

As expected, the main feature of tropospheric 0 3 during this 
period is the widespread maximum (total column 0 3 40 DU or 
more) extending from southern equatorial Africa across the 
Atlantic to South America. Examination of averaged winds 
and satellite-derived fire locations strengthens the case for a 
strong biomass burning component of this 0 3 feature. North of 
the equator regions of elevated tropospheric 0 3 are less uni- 
form, which is consistent with biogenic sources, lightning at the 


ITCZ, and lower levels of biomass burning (rainy season). 
Elevated 0 3 levels near Australia and New Guinea may signal 
biomass burning in that region. 
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